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Abstract 

We present a computation of the cross section for inclusive Higgs production in gluon-gluon 
fusion for finite values of the top mass in perturbative QCD to all orders in the limit of high 
partonic center-of-mass energy. We show that at NLO the high energy contribution accounts 
for most of the difference between the result found with finite top mass and that obtained in 
the limit mt — > 00. We use our result to improve the known NNLO order result obtained at 
rut — 00. We estimate the effect of the high energy NNLO mt dependence on the K factor to 
be of the order of a few per cent. 
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1 The cross section in the soft limit and in the hard limit 



The determination of liiglier-order corrections to collider processes [1], and specifically Higgs 
production [2] in perturbative QCD is becoming increasingly important in view of forthcoming 
phenomenology at the LHC The dominant Higgs production mechanism in the standard model 
is inclusive gluon-gluon fusion {gg —>■ H + X) through a top loop. The next-to-leading order 
corrections to this process were computed several years ago [3,4] and turn out to be very large 
(of order 100%). The bulk of this large correction comes from the radiation of soft and coUinear 
gluons [5], which give the leading contribution in the soft limit in which the partonic center-of- 
mass energy s tends to the Higgs mass m'jj, and which at LHC energies turns out to dominate 
the hadronic cross section after convolution with the parton distributions. 

This dominant contribution does not resolve the effective gluon-gluon-higgs (ggH) coupling 
induced by the top loop. As a consequence, the NLO correction can be calculated [6, 7] quite 
accurately in the limit — oo, where it simplifies considerably because the ggH coupling 
becomes pointlike and the corresponding Feynman diagrams have one less loop. Recently, the 
NNLO corrections to this process have been computed in the — > oo limit [8]. The NNLO 
result appears to be perturbatively quite stable, and this stability is confirmed upon inclusion [9] 
of terms in the next few orders which are logarithmically enhanced as s — mn, which can be 
determined [10] using soft resummation methods. This suggests that also at NNLO the large 
mt approximation should provide a good approximation to the yet unknown exact result. 

However, this is only true for the total inclusive cross section: for example, if one looks at the 
production of Higgs plus jets, if the transverse momentum is large the infinite approximation 
fails [11]. Indeed, even though the mf-independent contribution from soft and coUinear radiation 
turns out to dominate the cross section at the hadronic level, it does not necessarily provide a 
good approximation to the partonic cross section in a fixed kinematical region. In particular, 
the infinite approximation, which becomes exact in the soft limit, fails in the opposite (hard) 
limit of large center-of-mass energy. This is due to the fact that the ggH vertex is pointlike 
in the infinite rrit limit, whereas for finite rrit the quark loop provides a form factor (as we 
shall see explicitly below). Clearly, a point-like interaction has a completely different high 
energy behaviour than a resolved interaction which is softened by a form factor: in fact one can 
show [12] that a point-like interaction at n-th perturbative order has double energy logs while 
a resolved interaction has only single logs. 

This means that as s — oo the gg H + X partonic cross section a behaves as 



Hence, as the center-of-mass energy grows, eventually rrit ^ oo ceases to be a good approxi- 
mation to the exact result. It is clear from eq. ([1]) that this high energy deviation between the 
exact and approximate behaviour is stronger at higher orders, so one might expect the relative 
accuracy of the infinite rut approximation to the A;-th order perturbative contribution to the 
cross section to become worse as the perturbative order increases. Conversely, this suggests that 
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it might be worth determining the high energy behaviour of the exact cross section and use the 
result to improve the infinite rrit result, which is much less difficult to determine. Eventually, a 
full resummation of these contributions might also become relevant. 

The leading high energy contributions to this process in the infinite rrit limit have in fact 
been computed some time ago in Ref. [13]: this amounts to a determination of the coefficient of 
the double logs eq. ([I]), in the pointlike case. In this paper, we compute the coefficients of the 
single logs eq. ([1]) in the resolved (exact) case. Our result takes the form of a double integral, 
whose numerical evaluation order by order in a Taylor expansion gives the coefficient of the logs 
eq. ([I]) (at the lowest perturbative order the integral can be computed in closed form). After 
checking our result against the known full NLO result of Ref. [3,4], we will discuss the way 
knowledge of the exact high energy behaviour of the cross section at a given order can be used 
to improve the infinite rrit result, using the NLO case, where everything is known, as a testing 
ground. We will show that in fact, at NLO the different high energy behaviour eq. ([1]) accounts 
for most of the difference between the exact and infinite rrit cross sections. We will then repeat 
this analysis in the NNLO case, where only the infinite rrit result is currently known. We will 
show that in fact at this order the contribution of the logarithmically enhanced terms which 
dominate the partonic cross section at high energy is substantial even for moderate values of the 
partonic center-of-mass energy, such as s ~ 2m^. 

The calculation of the leading high energy logs is presented in section 2, while in section 3 
we discuss its use to improve the NLO and NNLO results. The appendix collects the explicit 
expressions of the form factors which parametrize the amplitude for the process gg H with 
two off-shell gluons, which is required for the calculation of sect. 2. 

2 Determination of the leading high energy logarithms 

2.1 Definitions, kinematics and computational procedure 

We compute the total inclusive partonic cross section a{gg H + X) in an expansion in power 
of a^, as a function of the partonic center-of-mass energy s: 



where the dimensionless variables r and yt parametrize respectively the partonic center-of-mass 
energy and the dependence on the top mass: 
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The corresponding contribution to the hadronic cross section a can be obtained by convolution 
with the gluon-gluon parton luminosity C: 
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where ghi{xi, Q"^) is the gluon distribution in the i-th incoming hadron and in eq. 
sionless variables Th parametrizes the hadronic center-of-mass energy s 



(5) 
(6) 

the dimen- 

(7) 



Note that < r/i < r < 1, and that ii yt < \ then the intermediate tt pair produced by the 
gluon-gluon fusion can go on shell. 

It is convenient to define a dimensionless hard coefficient function C(as(m|^); r, y^) 



gg {as;T;yt,mjj) 
C(a,(m^),r, yt) 



(To(i/t)C(as(m^),r, yt) 
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where (ToS{1 — t) is the leading order cross section, determined long ago in ref. [14]: 
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We also define the Mellin transform 
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C{as{mj,),N,yt)= / dr T^-'C{as{mj,),T,yt) 
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denoted with the same symbol by slight abuse of notation. 

We are interested in the determination of the leading high energy contributions to the partonic 
cross section a{gg H + X), namely, the leading contributions to C{as{'m'jj),T,yt) as r — 
to all orders in asirrijj). Order by order in ^^(m^), these correspond to the highest rightmost 
pole in in the expansion in powers of as(m|^) of C{as{'rn\j),N,yt). The leading singular 
contributions to the partonic cross section (T{gg ^ H + X) to all orders can be extracted [12] 
from the computation of the cross section for a slightly different process, namely, the cross 
section cros{gg — > H) computed at leading order, but with incoming off-shell gluons, a suitable 
choice of kinematics and a suitable prescription for the sum over polarizations. 

The procedure used for this determination is based on the so-called high energy (or kt) 
factorization [12], and consists of the following steps. 
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One computes the matrix element A^(^^(/ci, ^2) for the leading-order process gg ^ H at 
leading order with two incoming off-shell gluons with polarization indices n, v and color 
indices a, h. The momenta /ci, ^2 of the gluons in the center-of-mass frame of the hadronic 
collision admit the Sudakov decomposition at high energy 

ki = ZiPi + ki, (13) 

where Pi are lightlike vectors such that Pi ■p2 7^ 0, and kj are transverse vectors, kj ■ Pj = 
for all The gluons have virtualities 

fcNk? = -|k.r. (14) 

The cross section Cosigg ~^ H) is computed averaging over incoming and summing over 
outgoing spin and color: 

^ofT = (^1) E^''(^2)e*^?(A:2) ^p, (is) 

Ai A2 

where the flux factor 

J = 2{h ■ A;2 - ki ■ k2) (16) 
is determined on the surface orthogonal to pi, p2 eq. ( fT3l) . and the phase space is 



^'^-^^fi-l-^l. (IT) 



Note that the kinematics for a 2 ^ 1 process is fixed, so eq. f|T5l) gives the total cross 
section and no phase-space integration is needed. 

The sums over gluon polarizations are given by 

Here, the virtualities will be parametrized through the dimensionless variables 

= (19) 



The reduced cross section a, obtained extracting an overall factor mjj, 

mjjaosigg ^ H) = a{yt; ^1, 6, V^, ^), (20) 

is then a dimensionless function a{yt;^i,^2,^,z) of the parameter yt eq. (jlj) and of the 
kinematic variables ^1, ^2, the relative angle <f of the two transverse momenta 

^ = cos- ( ^] , (21) 



and 



|ki||k2 



my my 



2ziZ2Pi ■ P2 2 {ki ■ A;2 - ki • k2) 
Note that, in the collinear limit ki, k2 ^ 0, 2; eq. (l22l) reduces to r eq. 
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The reduced cross section is averaged over and its dependence on z eq. ( 1221) is Mellin- 
transformed: ^ ^ 

a(iV,ei,6)= / rf^^"^-' [^^a{yu^i,^2,^,z). (23) 



27r 
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The dependence on is also MeUin-transformed, and the coefficient of the coUinear pole 
in Ml, M2 is extracted: 



00 



h{N,M^,M2) = M^M2 d^i d^2er~'er HN,^i,^2). (24) 

Jo Jo 

Note that the integral in eq. ( l24l) has a simple pole in both Mi = and M2 = 0. The residue 
of this pole is the usual hard coefficient function as determined in coUinear factorization, 
which is thus C{N) = h{N, 0, 0). 



The leading singularities of the hard coefficient function eq. (1121) are obtained by expanding 
in powers of as at fixed ag/N the function obtained when Mi and M2 in eq. are 
identified with the leading singularities of the largest eigenvalue of the singlet anomalous 
dimension matrix, namely 

mlao{yt)Ciasiml),N,yt) = h (iV,7, (^) ,7, (^)) [1 + («,)] . (25) 

Here, 7^ is the leading order term in the expansion of the large eigenvalue 7^ of the singlet 
anomalous dimension matrix in powers of at fixed ag/N: 

7+KiV)=7s(^)+7..(|)+..., (26) 



with [15] 



\n)~ 2^'''' \ nN J 

n=l ^ / 



^s^)=ycn ^] ; c. = l,0,0,2C(3),..., (27) 



where Ca = 3. 

So far, this procedure has been used to determine the leading nontrivial singularities to the 
hard coefficients for a small number of processes: heavy quark photo- and electro-production [12] , 
deep-inelastic scattering [16], heavy quark hadropro duct ion [17,18], and Higgs production in the 
infinite rrit limit [13]. 

2.2 Cross section for Higgs production from two off-shell gluons 

The leading-order amplitude for the production of a Higgs in the fusion of two off-shell gluons 
with momenta ki and ^2 and color a,b is given by the single triangle diagram, and it is equal to 
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where the strong couphng is = |^ and the top Yukawa couphng is given hy ht = ^ in 
terms of the Higgs vacuum-expectation value v, related to the Fermi coupling by Gp = 
The dimensionless form factors Ai{^i,^2]yt) and ^2(^1; ^2; ?/t) have been computed in ref. [11]; 
their explicit expression is given in the appendix. They were subsequently rederived in Ref. [19], 
where an expression for the Higgs production cross section from the fusion of two off-shell gluons 
was also determined, but was not used to obtain the high energy corrections to perturbative 
coefficient functions. 

The spin- and colour-averaged reduced cross section eq. (l20l) is then found using eq. (fTSll . 
with the phase space eq. (fTTI) . We get 



2z 



5 Q - 1 - ^1 - ^2 - V6i2COSV9^ . 



(29) 



Because of the momentum-conserving delta, the Mellin transform with respect to z is trivial, 
and the reduced cross section eq. ( l23i) is given by 



+ 6)^(l + v^cos¥.) 



N 



(30) 



|Ai|2cOsV + 661^31') + 



where we have defined the dimensionless variable 

466 



+ 6 + 6) - (^!^2 + AiA;)] cos if 



a 



:i+ei + 6)^' 



(31) 



The three form factors Ai are independent of if, so all the angular integrals can be performed 
in terms of hypergeometric functions, with the result 



1^1 



;i+6 + 6 



^ ^ ,N N + 1 
2Fi(-,^— ,2,a) 
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N[\A^\^il + ^,+^2)-iAlA, + A^A; 
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A^ + 1 



1, a] 



1 



1 + 6 + 6 



2 2 

N + 1 N + 2 



2, a)}. (32) 



In the limit 00, using the behaviour of the form factors eq. (!60|) the term in square 

brackets in eq. (132!) as well as the term proportional to ^3 are seen to vanish. The remaining 
terms, proportional to Ai, give the result in the pointlike limit. The reduced cross section in 
this limit was already derived in ref. [13] (see eq. (9) of that reference): our result differs from 
that of ref. [13], though the disagreement is by terms of relative 0{N), hence it is immaterial 
for the subsequent determination of the leading singularities of the hard coefficient function. 
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2.3 High energy behaviour 

The leading singularities of the coefficient function can now be determined from the Mellin 
transform h{N, Mi, M2) eq. ( !24l) of the reduced cross section eq. ( !32l) . letting Mi = M2 = 
7<;(a<i/iV) according to eq. ( l25l) . and expanding in powers of as (i.e. effectively in powers of 
Ml, M2) and then in powers of N about N = 0. In the pointlike case {rrit 00) the Mellin 
integral eq. (!24|) diverges for all Mi, M2 when = 0, and it only has a region of convergence 
when > 0. As a consequence, the function h{N, Mi, M2) eq. (12^ has singularities in the M 
plane whose location depends on the value of A^, namely, simple poles of the form j^ztm^ztm^'- the 
expansion in powers of Mj has finite radius of convergence Mj < A^, leading to an expansion in 
powers of ^ and thus double poles when Mj = 7^. 

In the resolved case (finite rrit) we expect the Mellin integral to converge when A^ = at least 
for < Mi < Mq, for some real positive Mq. We can then set A^ = 0, and obtain the leading 
singularities of the coefficient function from the expansion in powers of M of h{0, M, M), letting 
M = 7s. This turns out to be indeed the case: when A^ = 0, ct(A^, ^1,^2) only depends on ^1, ^2 
through the form factors, and the combination of form factors which appear in a eq. f l32l) is 
regular when ^1, ^2—^0 (see eq. (l63|l ). while it vanishes when ^1, ^2 — ^ 00 (see eq. (|6^ ). Hence, 
we can let A^ = in a, and get 



h{0,Mi,M2) = %^-K^alGFm]jyl 



2 



p+oc r+oo 

XM1M2 / rfeief"' / d^2^^''-' 

Jo Jo 



(33) 



Because the term in square brackets in eq. (l33l) tends to a constant as ^1, ^2 — > 0, the 
integrals in eq. (l33ll have an isolated simple pole in Mi and M2, and thus the Taylor expansion 
of h{N, Ml, M2) has a finite radius of convergence. We can then determine the Taylor coefficients 
by expanding the integrand of eq. fl33l) and integrating term by term. It follows from eqs. fl25tl27p 
that knowledge of the coefficients up to A;-th order in both Mi and M2 is necessary and sufficient 
to determine the leading singularity of the coefficient function up to order a^. 

Let us now determine the leading singularities of ffist three coefficients of the expansion of 
the coefficient function eq. The constant term determines the leading-order result ctq eq. (jS]): 

mj,ao{yt) = h{0,0,0). (34) 

Using the on-shell limit of the form factors (see eq. (1^ of the appendix) in eq. we reproduce 
the well-known result eq. (fTOl) . 

The next-to-leading order term C^^\N,yt) is determined by noting that 

p+00 

/i(0,M,0) = aV27c' a'^Gprnly^M da'"" 0)|' 

Jo 

= h{0, 0, 0) - 8V27r'' a^G pmly^ M d^\n^^^^^^^^ + 0{M^). (35) 

Jo dt, 
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270 


1.1409 


4.1738 


280 


0.9568 


3.9828 


290 


0.7731 


3.8268 


300 


0.5884 


3.7049 



Table 1: Values of the coefficients eq. (l36ll and eq. (138|) of the 0{as/N) and 0((q;s/A^)^) of the 
leading singularities of the coefficient function C^asimj^); N,yt) eq. (fT2!) . 



Equations fl25ti27p then immediately imply that 



C«(iV,y,)=C«(y*)^[l + 0(iV)], 

C(^) ^ , r^eine^^#^. (36) 



(i-i(i-4y,).o(i/*)')r^o 

The value of the coefficient C^^\ determined from a numerical evaluation of the integral in 
eq. fl36l) . is tabulated in table 1 as a function of the Higgs mass. Upon inverse Mellin transfor- 
mation, one finds that 

limC«(r,yi) = C^CW(2/t)- (37) 

T— >0 

The values given in table 1 are indeed found to be in perfect agreement with a numerical 
evaluation of the small r limit of the full NLO coefficient function C^^\T,yt) [4], for which we 
have used the form given in ref. [20]. 

Turning finally to the determination of the hitherto unknown NNLO leading singularity, we 
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evaluate the O(M^) terms in the expansion eq. (^31) : by using again eqs. fl25H27l) we find 



C(2)(iV,y,)=C(2)(y,)^[l + 0(iV)] 



(87T 



2\2 



d^2 



In ^1 In^- 



+ 21^3(^1,6)1' 



(3^ 



The value of the NNLO coefficient C^'^\yt) obtained from numerical evaluation of the integrals 
in eq. (138!) is also tabulated in table 1. This is the main result of the present paper. 



3 Improvement of the NLO and NNLO cross sections 

Knowledge of the leading small r behaviour of the exact coefficient function C{as{m'jj);T,yt) 
eq. ([8]) can be used to improve its determination. Indeed, as discussed in section 1, we expect 
the pointlike (m^ — oo) approximation to be quite accurate at large r, whereas we know that 
it must break down as r ^ 0. Specifically, the small r behaviour of the coefficient function is 
dominated by the highest rightmost singularity in C^asirrijj); N, yt) eq. f|T2l) . which for the exact 
result is a fc-th order pole but becomes a 2A;-th order pole in the pointlike approximation. Hence 
the pointlike approximation displays a spurious stronger growth eq. ([1]) at small enough r. 

Having determined the exact small r behaviour up to NNLO, we can improve the approximate 
pointlike determination of the coefficient function by subtracting its spurious small r growth and 
replacing it with the exact behaviour. We discuss first the NLO case, where the full exact result 
is known, and then turn to the NNLO where only the rrit oo result is available. 



3.1 NLO results 

At NLO the small r behaviour of the coefficient function in the pointlike approximation is 
dominated by a double pole, whereas it is given by the simple pole eq. fl36|) in the exact case. 
This corresponds to an exact NLO contribution C^^\T,yt) which tends to a constant at small r, 
whereas the pointlike approximation to it grows as Inr: 

C«(r,oo) = d'^Ur) + O (r) ; d'^^i^) = c\ln r + c\ (39) 
C«(r,y,) = rf«(r,y,)+0(r); rf«(r, y,) = 3C«(y,), (40) 

where C^^\yt) is tabulated in table 1, while from Refs. [4,6,7] we get 

c\ = -Q; c\ = -^. (41) 

The NLO term C^^\t, yt) eq. (IE]) is plotted as a function of r in fig. [T], both in the pointlike 
{rrit — > oo) approximation, and in its exact form computed with increasing values of the Higgs 
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Figure 1: The hard coefficient C^^\T,yt) eq. ([9]) (parton-level coefficient function normahzed to 
the Born result) plotted as a function of r. The curves from top to bottom on the left correspond 
to mj = oo (black), and to = 170.9 GeV (red), with = 130, 180, 230, 280 GeV. 



mass, i.e. decreasing values of yt- It is apparent that the pointlike approximation is very accurate, 
up to the point where the spurious logarithmic growth eq. ( |39ll sets in. 

We can construct an approximation to C^^\T,yt) by combining the pointlike approximation 
with the exact small r behaviour: 

C(^)'^PP-(r,y,) ^ C(^)(r,oo) + [4^)(r,y,) - d^;U^)] T{t) (42) 

where d'^{T,yt) and d^^^^^{T) are defined as in eq. pUj) and eq. fl39|) respectively, while T(r) is 
an interpolating function, which we may introduce in order to tune the point where the small r 
behaviour given by rfcx (t, yt) sets in. Clearly, as r ^ the approximation eq. fH2l) reproduces the 



exact small r behaviour of the exact coefficient function eq. ( HOj) . provided only the interpolating 
function limT-^o T{t) = 1. Furthermore, as discussed in section 1, the behaviour of the coefficient 
function C^^\T,yt) as r 1 is to all orders controlled by soft gluon radiation, which leads to 
contributions to C^^\t, yt) which do not depend on yt and diverge as r — > 1. Hence, the pointlike 
approximation is exact as r ^ 1. Because the functions d'^{T,yt) and d^^^^^{T) are regular as 



r — * 1, this exact behaviour is also reproduced by the approximation eq. (142]) . provided only 
lim^^i T(r) is finite. Hence, C'(^)''^PP-(r, yt) reproduces the exact C^^\t, yt) as r ^ up to terms 
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Figure 2: The hard coefficient C(i)(r,?/j) eq. with = 130 GeV (left) and = 280 GeV 
(right). The sohd curves correspond to = oo (black) and = 170.9 GeV (red), (same 
as fig. [1]). The three blue curves correspond to the approximation eq. (142 f 43 1) , with k = 
(dot-dashed), k = 5 (dotted), k = 20 (dashed). 



that vanish as r — and as r — > 1 up to terms that are nonsingular as r — 1, even when 
r(r) = 1. 

Nevertheless, we may also choose T(r) in such a way that T{1) = (while T(0) = 1 always), 
so that C^^\T,yt) agrees with the pointlike approximation C^^\t, oo) in some neighborhood of 
r = 1. For instance, we can let 

T(r) = (1 - r)^ (43) 

with k real and positive, so that the first k orders of the Taylor expansion about r = 1 of 
C^^'''^PP (r, i/i) and the pointlike approximation coincide. By varying the value of fc, we can 
choose the matching point tq, such that C^^^'^^'^'(t, yt) only differs significantly from the pointlike 
approximation if r < tq: a larger value of k leads to a smaller value of tq. 

In fig. 2 we compare the approximate NLO term eq. fH2]) to the exact and pointlike results, 
for two different values of yt, with T(r) given by eq. fH51) and a choice of k which leads to different 
values of the matching between approximate and pointlike curves. It appears that an optimal 
matching is obtained by choosing k in such a way that the approximation eq. fH2l) matches 
the pointlike result close to the point where the logarithmic growth of the latter intersects 
the asymptotic constant value of the exact result. Note that this optimal matching could be 
determined without knowledge of the exact result. With this choice, the approximation eq. fH2]) 
differs from the exact result for the NLO contribution to the partonic cross section by less than 
5% for all values of r. 
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Figure 3: The hard coefficient C^'^\T,yt) eq. ([9]) (parton-level coefficient function normahzed 
to the Born result) plotted as a function of r. The curves from top to bottom on the left 
correspond to = oo (black), and to the approximation eq. fHTl) with T(r) eq. fH3|) and k = 5, 
and mt = 170.9 GeV (red), with tuh = 130, 180, 230, 280 GeV. 



3.2 NNLO and beyond 

At NNLO, the pointlike approximation to the coefficient function has a quadruple pole at = 0, 
corresponding to a In^ r rise, while the exact result only has a double pole, and thus it rises 
linearly with Inr: 

C(2)(r,oo) = 41,(r) + 0(r°); S^ti^) = c\\n' r + c\\n' r + cMnr (44) 
C^'\r,y,) = 4^)(r,y,)+0(r°); S^ir^y,) = -9C^'\y,)\nT, (45) 

where C^'^\yt) is tabulated in table 1, while from Ref. [8] we get 

2 „ 2 231 17 2 f 2333 „ o\ 641 



c^4 = -6; c 3 = +nf — - c 2 = + 37r^ +nf , (46) 

•3 4 ^18' V 8 / 108' ^ ' 

where rij the number of flavors. 

At this order, the exact form of C^'^\T,yt) is not known. However, analogously to the NLO 
case, we construct an approximation to it based on its determination [8] in the pointlike limit. 



12 



1500 



1000 



500 



-1 — I I I 1 1 1 1 1 — I — I I I 1 1 1 1 1 — I — I I I 1 1 1 1 1 — I — I I I 1 1 1 







5 CO 



I 



c(^)(t) 



200 



100 — 



_i I III' 



10" 



10" 



10" 



10" 



10" 




Figure 4: The hard coefficient C^'^\T,yt) eq. ([9]) with uih = 130 GeV, plotted versus r on a 
logarithmic (left) or linear (right) scale. The solid black curve corresponds to rrit = oo (black, 
same as fig. [3])), and the three blue curves are the approximation eq. ( H7|) with rrit = 170.9 GeV 
and T(r) eq. (143|1 with and k = (dot-dashed), k = 5 (dotted, same as fig. [3]), k = 20 (dashed). 



combined with the exact small r behaviour eq. ( l38i) : 

C(^)'^PP-(r, y,) ^ C^'\t, oo) + [4^J(r, y,) - <l,(r)] T(r) (47) 

fr)\ fr)\ _____ 

with dcx{T,yt) and rfpoint(''") defined in eq. (H5ll and eq. (I44l) respectively, and T(r) an inter- 
polating function as discussed in section 3.1. Note that as r — the approximation eq. ( HTll 
only reproduces the exact result up to a constant, whereas at NLO the approximation eq. fH2|) 
reproduces the exact result up to terms which vanish at least as 0(r) . 

The approximation to the exact result C^'^^'^'^^'{T,yt), computed using C*-^^ from table 1 with 
four different values of the Higgs mass, and taking T(r) eq. ( H3l) with A; = 5 is compared in fig. 3 
to the pointlike approximation C^'^^'^^^'{T,yt) of ref. [8] (with Uf = 5). In figures 4-5 we further 
compare the results obtained with different choices of the matching function T(r) eq. (H3ll . and 
the same two values of the Higgs mass used to produce figs. 2-3 at NLO. 

At this order, the contribution from the leading small r logs to the pointlike C^'^^'^^^'{t, oo) 
is sizable even for large r. Indeed, figs. 4-5 show that the behaviour of C*-^^ around its local 
maximum at r ~ 0.65 receives a sizable contribution from the Inr rise and In^ r drop eq. (1441) . 
If these are removed by using eqs. fli71ll3D with = 0, the shape of C^^^ around the maximum is 
affected significantly, but if the matching is moved to smaller r by choosing A; ~ 5 the maximum 
is reproduced. Hence, whereas we can still obtain a rather smooth matching at any desired value 
of T the choice of the optimal value of r is not obvious. In particular, matching at a value of r 
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Figure 5: Same as fig. HJ but now with rriH = 280 GeV. 



where the contribution of the asymptotically spurious In^ r becomes significant leads to rather 
large values of the matching point r ~ 0.6. Anyway, it is clear that the pointlike approximation 
breaks down for r ^ 0.1. 

Contributions beyond NNLO in the expansion of h{N,'-fs,'js) eq. fl25l) in powers of ^ can be 
determined by pursuing the expansion of h{0, M, M) eq. fl33l) in powers of M, and determining 
numerically the ensuing integrals, which have the form of eqs. ( !36]l38l) . but with higher order 
powers of In^i, In ^2- The series of contributions to the coefficient function eq. (IHl) thus obtained 
has a finite radius of convergence in A^-Mellin space, dictated by the location of the rightmost 
singularity in 7^, and thus in r space it converges for all < r < 1 [22]. Therefore, its resum- 
mation can be accomplished to arbitrary accuracy by computation of a finite number of terms. 
This resummation, however, induces spurious singularities in the A^-space coefficient function, 
which can be removed by the inclusion of a suitable class of formally subleading running-coupling 
corrections, as recently shown in Ref. [23]. 

3.3 K factors 

The accuracy of the various approximations at the level of hadronic observables clearly depends 
on the individual process. For the total inclusive cross section eq. (I5l), as is well known, the 
pointlike approximation is actually very good, and thus the impact of the improvement eq. P2|) 
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^NLO 


^NNLO 


rriH — 


130 GeV 




pointlike 




36.69 




OOO 


exact 




36.58 




n.a. 


appr., k = 5 




37.64 






appi . , ft — zu 




36.66 




uoo 


rriH = 


r 


80 GeV 




pointlike 




38.08 




716 


exact 




37.47 




n.a. 


appr., k = 2 




37.97 




670 


appr., k = 5 




37.73 




693 



Table 2: The NLO and NNLO contributions to the K factor eq. (HHj) . computed with center-of- 
mass energy s = 14 TeV, and nit oo, denoted with pointlike, or = 170.9 GeV, denoted 
with exact or approximate. The approximate result uses eqs. fl42|47l) . with T(r) eq. and 
the value of k given in the table. The MRST2002 [21] gluon distribution has been used. 



is moderate. To give a quantitative assessment, we define a K factor by letting: 

+ OK(m|)), (48) 

where a^^ is the leading-order form of the contribution eq. ([5]) of the gluon-gluon channel to 
total hadronic cross section. The value of the NLO and NNLO contributions to the K factors, 
determined using the MRST2002 [21] gluon distribution in eq. ([5]) are given in table 2 at LHC 
energies for two values of the Higgs mass, both in the pointlike, exact and approximate (eq. ( H2l) 
and eq. (jTTl) ) cases. 

At NLO with rriH = 130 GeV ("light"), the pointlike approximation to k^^^ deviates by 
0.3% from the exact result, and even with rnn = 280 GeV ("heavy") it only deviates by 1.6%. 
It should be kept in mind, however, that itself is quite large: for ~ 0.1, it amounts 

to a ~ 100% contribution to the K factor eq. (jUD- Hence, the error made using the pointlike 
approximation is between the per mille and the per cent level, and thus not entirely negligible 
in a precision analysis. 

Using the approximation eqs. (142ti43p with the values k = 20 for light Higgs and k = 5 for 
heavy Higgs, which are seen from fig. 2 to give good matching, the deviation can be reduced to 
0.2% and 0.7% respectively, and even more accurate results could be obtained by an optimization 
of the matching. However, a poor choice of the matching (such as /c = 5 for light Higgs or k = 2 
for heavy Higgs) can lead to a result at the hadronic level which is actually closer to the pointlike 
approximation, or even worse than it. It is clear that at the partonic level the small r behaviour 
eq. fl39l) accounts for most of the discrepancy between the exact and pointlike results, and even 
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the determination of a hadronic observable which depends very httle on the parton-level small 
r behaviour can be improved very substantially for values of relevant for LHC by using the 
approximation eq. fH2l) . 

The NNLO contribution not known. Its values computed in the pointlike approxi- 

mation, or with the approximation eqs. ( I47f43i) and different choices of k are shown in table 2. 
Even at the inclusive hadronic level, now the size of the NNLO contribution can change up 
to about 5 — 10% if the matching is performed at large r. Furthermore, also quite 

large: with ^ 0.1, it amounts to a ~ 50% correction to the leading order, and thus to a 
further ~ 25% correction to the K factor. Therefore, the impact of the pointlike approximation 
at NNLO is up to several per cent of the total K factor, rather larger that the impact of the 
pointlike approximation at NLO, and comparable to uncertainties which are currently discussed 
in precision studies at NNLO. 

4 Outlook 

2 

In this paper we have determined the leading high energy (i.e. small r = singularities of the 
cross section for Higgs production in gluon-gluon fusion to all orders in the strong coupling, by 
providing an expression (eq. fl33l) ) whence the coefficients of these singularities can be obtained by 
Taylor expanding and computing a double integral. We have given explicit numerical expressions 
for these coefficients up to NNLO. 

The high energy behaviour of this cross section is different according to whether it is deter- 
mined with finite mt or with — ^ oo (pointlike approximation). It turns out that at NLO this 
different high energy behaviour is responsible for most of the discrepancy between the pointlike 
approximation and the exact result. As a consequence, an accurate approximation to the exact 
result can be constructed by combining the pointlike approximation at large r with the exact 
small r behaviour. Some care must be taken in matching, but very accurate results can be 
obtained by simply choosing the matching point as that where the spurious small r behaviour 
of the pointlike behaviour sets in. 

At NNLO, where the exact result is not known, the impact of the high energy behaviour 
turns out to be large even for moderate values of r ~ 0.5. Hence, an approximation constructed 
analogously to that which is successful at NLO, namely matching the pointlike limit to the 
asymptotic exact behaviour at the point where the asymptotically spurious terms become signif- 
icant, leads to an approximate result which differs significantly from the pointlike approximation 
for most values of the partonic center-of-mass energy. 

The effect of these high energy terms on the total inclusive hadronic cross section remains 
quite small, because the latter is dominated by the region of low partonic center-of-mass energy, 
partly due to shape of the gluon parton distributions, which are peaked in the region where the 
gluons carry a small fraction of the incoming nucleon's energy, and partly because the partonic 
cross section is peaked in the threshold r ^ 1 region. Even so, the pointlike determination of 
the NNLO contribution to the total hadronic cross section can be off by almost 5-10% due to 
this spurious high energy behaviour, especially for relatively large values of tuh ~ 200 GeV. 
Because the NLO and NNLO corrections to the cross section are quite large, the overall effect of 
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these terms on the cross section is at the per cent level, and in particular their effect at NNLO 
is rather larger than at NLO. 

A study of the phenomenological implications of these results is thus relevant for a precision 
determination of the Higgs production cross section. 
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A Form factors 



The form factors in eq. (!28l) are given by 



+ 



A3 

[i?o(-6)-i?o(l)] 

[i?o(-6)-5o(l)] 
2 1 



A3 



1 + 6 + 6) -1-^ + 12^(1 + 6 + 6) 



26 , ,.66 
A3 -3 
26 , ,,66 
A3 



12^(1 + 6-6) 



+ 12:^(1-^1 + 6) 



;i + 6 + 6j 



A3 (47r)2 

^2(6,6,1/*) = Co(6,6,2/0 

+ [i?o(-6)-i^o(l)] 
+ [i?o(-6) - fio(l)] 



2i/t-^(i + 6 + 6) + ^ 

^(1-6 + 6) 
^3 



1^(1 + 6-6 
^3 



+ 



with 



A3 = 1 + 6' + 6' - 266 + 2(6 + 6) = (1 + 6 + 6)' - 466 • 

It is also convenient to define the form factor 



^3(6,6,2/t) 

The scalar integrals Bq and Cq are 

Boip) = 



66 



1+6 + 6 



A1-A2 



1 A yt - p 

tan 



P 



Boip) 



87r2 

1 Ip - Ayt 



p 



4:yt- p 



, if < p < 4?/, 



t , 



IGtt^ V P 



In 



p-^yt 



if p < or p > 4?/, 



t , 



(49) 



(50) 
(51) 
(52) 



(53) 
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where 



+U2{y+St) + U2{y-St) - Li2(y+5i ) - U2{y-6^ ) 
+Li2(a:+(5^) + Li2(;r_5^) - Li2(a;+52") - Li2(x_52") 

+Li2(2;+5^) + Li2(2;_5;^) - U2{z+6^) - Li2(2;_53 ) | , 



. _ -6 +6-1 r _ - 6 - 1 . _ + 6 

Ol = ?= , 02 = 7= , Os 



and 



2yt V V e 



;2 



2» V V 6 



y± = - 

1 

In the infinite top mass hmit the scalar integrals become 



yt-*oo 

so that the form factors reduce to 



lim m^^Ai = mjj-^; lim 4m^yl2 = mjj ^\ ^ ^ 

mt-^oo 4o7r mt ^oo 4o7r Z 



These limits also imply that 



mt—» oo 

In the on-shell limit the scalar integrals are 

1 



lim m^Az — 0. 
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so that 



^i(0,0) 
^2(0,0) 



1 1 



^7r2 327r2 



In^ 



1 1 



+ 



167r2 327r 



(%t - 1) 



(63) 



The high energy hmit of the form factors is trivially determined when ,^1 ^ 00, ,^2 ^ 00 with 
6 7^ 6: 



lim Ai(^i,^2,yt) = 0; lim ^3(^1, ^2, = 0; lim ^2(6, 6, ^t) = 7^^• 



(64) 



If ^1 ^ 00, ^2 ^ cxD with ^1 = ,^2 the limit is more subtle. In this case we get 



lim Ai(e,e,i/t) 



lim 



^ -ln--l + v/4yt-ltan-\/ ^ 



where we have defined 



Co(ei,6,l/t) =C'o(el,6,Z/t)vA3• 



However, it turns out that 



lim Co(E,C,Vt) = 7r-F= 



21n - -4 + 4v/4t/t- Itan^i 

e ^ V%-1 



hence we conclude that eq. fl64j) holds also when ^1 = 



(65) 
(66) 



O ( I ) , (67) 
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